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Abstract 

Using the monomer-dimer representation of the lattice Schwinger model, with Nf = 1 
Wilson fermions in the strong-coupling regime {(3 = 0), we evaluate its partition function, Z, 
exactly on finite lattices. By studying the zeroes of Z{k) in the complex plane {Re{k), Im{k)) 
for a large number of small lattices, we find the zeroes closest to the real axis for infinite 
stripes in temporal direction and spatial extent S = 2 and 3. We find evidence for the 
existence of a critical value for the hopping parameter in the thermodynamic limit S -^ oo 
on the real axis at about kc — 0.39. By looking at the behaviour of quantities, such as the 
chiral condensate, the chiral susceptibility and the third derivative of Z with respect to 1/2/c, 
close to the critical point kc, we find some indications for a continuous phase transition. 



1 Introduction 

The Schwinger model [|l|, that is 2-diniensional QED (QED2) with massless electrons, has 
always attracted the interest of theoreticians not only because it is a rather simple model, 
which can be solved analytically, but also (and mainly!) because many of its properties are 
quite similar to those of 4-dimensional QCD. Already in the basic version with Nf = 1 
fermion flavours, one recovers a lot of QCD-like properties such as confinement for fermions, 
chiral symmetry breaking, due to the anomaly in the U{1) axial current, and charge screening 

Ban- 

By virtue of this similarity, one is also tempted to consider the lattice version of the 
Schwinger model as a test-model for lattice 4-dimensional QCD (QCD4). Also here one 
faces the problem of the choice of a lattice scheme for fermions: the two most common 
choices are the well-known Wilson fermions ^ and the staggered (or KogutSusskind [KS]) 
fermions [§, ^]. Most of the lattice calculations done up to now for the Schwinger model 
used the staggered fermion formulation [0, H, P] , in which case the chiral limit is obtained by 
simply setting the bare fermion mass parameter m appearing in the lattice Lagrangian to 
zero: m — > 0. All these lattice calculations seem to reproduce well the expected properties 
of the continuum massless Schwinger model, known from analytical results. 

On the contrary, very little is known about the lattice Schwinger model with Wilson 
fermions. Our interest in the Wilson formulation of QED2 comes from the recently made 



observation ||T0| that the critical point in the hopping parameter, at which the chiral limit is 
reached, may not agree with the naive expectation (see below, for a more detailed discussion 
of this point). This may also be of relevance for an understanding of the complicated phase 



diagram found in the Wilson formulation of lattice QCD4 with large number of flavours [|T^ 

In the lattice action with Wilson fermions the bare fermion mass m does not appear 
explicitly as in the lattice action with KS fermions. It contains as parameters the coupling 
constant (3 and the hopping parameter k, which is related to the bare fermion mass. Therefore 
in the Wilson fermion formulation there is the problem of defining a chiral limit, correspond- 
ing to m —>■ 0. In lattice QCD4 with Wilson fermions the chiral limit, at a given value of 
(3, is reached when the hopping parameter k approaches a certain critical value kc{l3), often 
defined as the value of k for which the pion mass M^ vanishes. In fact, in the chiral limit of 
QCD4 the pion becomes the Goldstone boson of the spontaneously-broken chiral symmetry 
and its mass is expected to vanish as ^/m when m — > 0. It tacitly is assumed that this 
defines a critical point which coincides with a critical point (zero) of the partition function. 
It is expected that the situation is similar for the lattice Schwinger model, in the sense 
that there will be a critical point for each value of the gauge coupling, /c = fcc(/3), which 
defines the chiral limit. The continuum chiral limit will be reached following this line up 
to kc{l3 = 00) = l/2d = 1/4. However, we cannot determine kc{(3) in the same way as 
in QCD4, since in the Schwinger model with Nf = 1 we have no Goldstone boson in the 
chiral limit m — *^ 0: the ^(1) chiral symmetry is broken by the anomaly and one is left with 
a massive pseudoscalar, similar to the rj' in QCD4. A determination of kc{l3) thus has to 
proceed through the direct investigation of the singularities of Z{k). 

A quite common attitude is to assume that kdP) coincides with the convergence radius 



k{l3) for the joint expansion in the hopping parameter k and the inverse gauge couphng 
f3. However, there is no proof that this is correct: for example, in the lattice Schwinger 
model with KS fermions we have that k < 1/2 [|T^, while kc = oo. For Wilson fermions, 
one does not know the precise values of fc(/5) and kc{(3). Our aim is to compare these two 
values at least in the strong coupling regime /3 = 0. It was already found in ref. |TD[ that 



^(0) < 1/2 and some indications from Monte Carlo simulations were reported, indicating 
that fcc(O) ^ A;(0). 

In this work we will determine kc{Q) by deriving analytically (using algebraic methods) 
the Lee- Yang zeroes [|12| of the partition function Z{k) for the lattice Schwinger model, with 



Nf = 1 Wilson fermions, in the strong-coupling regime (/5 = 0). For a finite lattice S x T 
these zeroes have a non-vanishing imaginary part, in the complex plane {Re{k),Im{k)), 
indicating that there is no real critical point for a finite lattice. This remains true for 
T —>■ oo with finite S. Yet, enlarging the lattice, they show a tendency to move towards the 
real fc-axis. By studying the zeroes in the complex plane of the partition function Z{k) for a 
large number of small lattices, and then extrapolating to the thermodynamic limit oo x oo, 
we will find evidence for the existence of a real critical value for the hopping parameter, at 
about kc — 0.39. We will also study some relevant quantities, such as the chiral condensate 
< ipip >, the chiral susceptibility and the third derivative, with respect to l/2k, of the 
partition function, in order to get some information about the question of the order of the 
phase transition. 

2 The method 

The action for the lattice Schwinger model is written as the sum of a gauge-action S'g[[/] 
and of a fermion-action Sf[iP,iI!,U]: 

S = SG[U] + SF[i;,ij,U]. (1) 

The gauge part 5'G[t/] is given by: 

SG[U]=(3j:[l-hup + Ul)], (2) 

p ^ 

where P = 1/e^, e being the usual electro-magnetic coupling constant. Since the gauge 
group is f/(l), the basic lattice gauge variable f/^(n), corresponding to the link connecting 
the sites n and n + fi, can be written in the form of a phase: 

Uf,{n) = exp[z0^(n)]. (3) 

In eq.(2) Up stands for the usual 1x1 Wilson plaquette, constructed using the link variables 

In the strong coupling limit (e^ —>■ oo), the coefficient (3 in front of ^^[f/] goes to zero 
and the total action 5* reduces simply to the fermion action Splipjip, U]. The action Sp for 



Wilson fermions (and only one flavour) can be written in the form: 

SF = ^Y.i'^^)Knm[U]i^{m), (4) 

^ n,m 

where k is the so-called hopping parameter, and each Knm-, for a given couple of lattice sites 
n and m, is a matrix in Dirac space: 

5nm ■ I -kJ2{i'^' - lt^)U^in)5n+f,,m + (r + 'J^)Ul{n - il)5n-f,,m]- (5) 

Therefore Knm is of the form: 

where the only non-vanishing matrices Mnm are those connecting neighbouring lattice sites: 

Mn,n+f.[U] = (r-7^)f/^(n), 

Mn,n-i.[U] = (r + 7^)f/;(n-/i). (7) 

The parameter r, which satisfies |r| < 1, is called the Wilson parameter. In the following we 
will consider only the case r = 1. The matrices 7^, with /i = 1, 2 are the 2x2 Euclidean Dirac 
matrices, corresponding to 1 + 1 space-time dimensions (in particular we will consider the 
index /x = 1 as corresponding to the time dimension and the index yU = 2 as corresponding 
to the space dimension): they satisfy the anticommutation relation {74,7^} = 26ikl2, with 
I2 being the 2x2 identity matrix. For our algebraic manipulations, we have chosen the 
following representation for the 7-matrices: 

-(J -°0 ■ -(? l)^ («) 

Finally, the partition function, in the strong coupling regime /3 = 0, is given by the following 
expression: 

Z{k) = j[D^D^] /[/)[/] e-^^[^'^'^l 



2k 



Y^i,{n)K^m[U]i^{ 



m] 



[Di)ij] / [DU]e "'"^ . (9) 

Following the standard normalization convention, we also eliminate the factor l/2k appearing 
in the exponent in eq.(9) by re-scaling the fermion fields with y2k: 



^ = V2k^, ^ = ^2k%lj. (10) 

When evaluating a matrix element of the form < n^=i i^A„{^n)i^Bniyn) >, in terms of the re- 
scaled fields, we must however keep in mind that the original correlation function is obtained 



by multiplying with (1/2/c)^. When considering a lattice with S lattice sites in the space 
direction and T lattice sites in the time direction, for a total oi N = S x T lattice sites, the 
partition function (9) becomes: 

Z{k) = (2A;)-2^Z(A;), (11) 



where we have defined: 



Z{k) = I [DipDip] I [DU]e "'™ . (12) 



It is exactly Z{k) that we have evaluated with algebraic methods for a large series of small 
lattices. First of all we have put Z{k) in a comfortable form for subsequent algebraic manip- 
ulations. Making use of the explicit expression (3) for the link variables U and remembering 
the Grassmann properties of the fermion fields, one finds the monomer-dimer representation 



1 1311 for the partition function Z{k): 

Z{k) = f[Di>DtP]l[F{n)L^{n)L2{n), (13) 

where F{n) is the monomer term at lattice site n, coming from the mass-term in the action 
{i.e. the bilinear diagonal term not containing the gauge fields U): 

F{n) = 1 - i)^{n)i)i{n) - tjj2{n)tjj2{n) + i!^{n)i!i{n)tjj2{n)tjj2{n). (14) 

(The indices 1 and 2 are the Dirac indices). The quantities L^{n)., with /i = 1,2, are the 
dimer terms defined on links of the lattice. They result from the direct integration over the 
gauge field U^{n) = exp[i(f)f^{n)], corresponding to the link n —> n + fi. Explicitly: 

J--K 2t[ 
where M^{n) and Nf^{n) are given by: 

Mf,{n) = k4j{n){l--f^)ij{n + ll), 

N^{n) = ki>{n + fi){l + ^f,)^{n). (16) 

Thanks to the particularly simple form (3) of the gauge variables U in the case of a U{1) 
gauge group, the integration in (15) can be performed in an elementary way, after having 
expanded as a power-series the first exponential. Making use of the explicit representation 
(8) for the 7-matrices and of the Grassmann properties of the fermion fields, one finds the 
following rather simple expressions for the one-link integrals Li{n) and L2{n): 

Li(n) = l + Mi{n)Ni{n) 

= l + 4:k'^ip2in)ip2in + i)^p^{n + i)ipi{n), 
L2{n) = l + M2{n)N2{n) 

= l + fc2^(n)(_^^ "^^)^(n + 2)-^(n + 2)Q J)^(n). (17) 



After inserting the expressions (14) and (17) in (13), one is left with an integral over the 
fermionic variables, which must be evaluated according to the integration rules for the 
Grassmann variables. Doing this "manually" turns out to be extremely boring and time- 
consuming, even for the very small 2x2 lattice. For this reason we have developed an 
algebraic method for evaluating the integral (13). We have used the algebraic computer 
language ^^Mathematica" and have implemented the basic rules for the Grassmann algebra: 



{Vi, Vj} = 0, / drii = 0, / dr]ir]k = 5, 



Hk- 



In this way we were able to write computer programs (for ^^Mathematica") for evaluating 
products of polynomials of arbitrary strings of Grassmann variables (like those appearing 
in eq.(13)) and for integrating them. In practice we have used the following strategy for 
calculating the partition function (13) for a given lattice of size S x T. First of all we have 
computed the transfer matrix, i.e. the products of all one-link-terms, Li{n) and L2{n), and 
all mass-terms F{n) belonging to a given space-like line (x, t) with x = 1,2, . . . , S. We call 
this object a ^^line" . It is a function of all Grassmann fields belonging to the sites of the line 

under consideration: 

s 

line[t, 0] = n ^i^, t)Li{x, t)L2{x, t). (19) 

x=l 

In evaluating this product we have already taken into account the torus-like topology of the 
lattice, in form of periodic boundary conditions for the Grassmann fields along space-like 
lines. The partition function is represented in terms of line\t, 0] as: 

_ T 

Z{k) = f[Di>Di>] f[ line[t, 0]. (20) 

-' t=i 

Starting with the object /me[t,0] one can now evaluate composite objects, such as the 
product of two adjacent lines: /me[t,0] ■/me[t + l, 0]. By virtue of eqs.(19) and (17), all other 
lines, different from line\t, 0] or line\t + 1, 0], do not depend on the Grassmann fields on the 
sites (x, t + 1), X = 1, 2, . . . , S*. Therefore we can integrate the product line\t, 0]-/me[t + l, 0] 
with respect to these Grassmann fields and obtain: 

line[t, 1] = / n #(^' ^ + ^)d'4'{x, t + l)]line[t, 0] ■ line[t + 1,0]. (21) 

•^ x=l 

We can then proceed in the same way and construct more extended objects. For example 
we can multiply line[t, 1] with line[t + 2,0], or even with line[t + 2, 1], and integrate over 
the Grassmann fields lying on the line {x,t + 2), x = 1,2, . . . ,S: In general we obtain 
line[t,l] by performing the integration over the (/) lines of intermediate Grassmann fields 
at X = t + 1, ...,t + /. Finally, for I = T — 1 the resulting object covers just the entire 
lattice S xT and we must again take into account the torus-like topology of the lattice, i. e. 
impose anti-periodic boundary conditions for the Grassmann fields along time-like lines (In 
fact, it turns out to be irrelevant, for the final result, if we impose periodic or anti-periodic 
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boundary conditions for these Grassmann fields). In the last step, we can integrate over the 
remaining Grassmann fields and obtain the final result for the partition function (13). 

It turns out that the CPU time required in this approach is entirely controlled by the 
initial spatial extent S of the lattice as this determines the number of Grassmann fields one 
can combine in a given string of fields. Doubling the length in time direction does not lead 
to a drastic increase of the computer time as the resulting object, line[t, 21] contains exactly 
the same number and types of strings of Grassmann fields as line[t,l]: the only additional 
complication results from the more complex structure of the coefficients of these strings, 
which become higher order polynomials in k for increasing T. As a curiosity, if we try to 
follow the strategy of ^^ extendinc}^ in the space-direction, instead of extending in the time- 
direction, we need (for our programs) CPU times which are orders of magnitude larger than 
in the previous case. This is simply due to the space-time asymmetry in the representation 
(8): but of course the partition function for a lattice S x T, namely Z{k,S,T), is exactly 
equal to the partition function Z{k,T,S) for a lattice T x S (One can always choose a 
representation for the 7 matrices in which 71 = 72 and 72 = 71). 

3 Results 

Following the computational strategy discussed in the previous section, we have evaluated 
the partition function Z[k), given by (13), for a large number of lattices of the form S x T 
with S* = 2 and 3 ranging from 2 x 2 up to 2 x 32 and from from 3 x 3 up to 3 x 16. These 
calculations could be performed on a workstation. For S* > 4 considerably more computer 
time and memory would be required. 

From the Grassmann properties of the fermion fields, it immediately follows that the 
function Z{k, S, T), for a given lattice having N = S x T sites, is a polynomial of order 2N 
in the hopping parameter k, 

N 

Z{k,S,T) = Y.a2nk'^, (22) 

n=0 

On a lattice with A^ lattice sites there are 4N different fermion fields (four fields for each 
site), so that, by virtue of the properties of the Grassmann algebra, one can construct strings 
of fermion fields with at most 4N fields. And since each power k'^ is always accompanied 
by four fermion fields, one can at most construct a string with 4A^ fields with a coefficient 
(proportional to ) /c^^ in front of it. 

We have computed the partition functions Z{k, S, T) for (S* = 2, T = 2, . . . , 32) and for 
(S* = 3,T = 3, . . . , 16): each of them is a polynomial of order 2A^ = 2(5* x T) in k, with 
ao = 1 {Z{k = 0,S,T) = 1) and a2n > for all n. The magnitude of these coefficients 
generally increases with the order of k, apart from deviations in the very first coefficients 
due the torus-like topology of the lattice which allows for special "trajectories" , made up 
of chains of links wrapping around the lattice. In Table 1 and Table 2 we report the list of 
these coefficients for the two lattices 2 x 32 and 3 x 16 respectively. Note that for S and 
T even the partition function Z{k,S,T) is a polynomial of the form Y,n=o (^ink^^ . This is 
due to the larger set of symmetries on such lattices, as will be discussed below. Let us first 



discuss the distribution of zeroes of the partition function. In [Fig.l] and in [Fig. 2] we show 
the distribution of the complex zeroes, {Re{k),Im{k)), of the partition function Z{k,S,T) 
for the 2 X 32 and 3 x 16 lattices, respectively. 

By virtue of (22), Z{k) is a polynomial in k'^ with real coefficients: so, if A; is a complex 
zero of Z{k), also —k and k* (the complex-conjugate of k) will be zeroes of Z{k). As a 
consequence of this, the distribution of zeroes {x,y) in the complex plane {Re{k),Im{k)) 
is invariant under the parity transformation {x,y) -^ {~x,—y) [P-symmetry) and under 
the comj'/ez-conju^'ate transformation {x,y) — * {x,—y) {C-symmetry) . In other words, the 
distribution of zeroes turns out to be symmetric under reflections with respect to the real 
and/or the imaginary /c-axis: this is evident from [Fig.l] and [Fig.2]. The distribution of 
zeroes in [Fig.l], for the lattice 2 x 32, has an additional symmetry under reflections with 
respect to the axis Re{k) —Im{k) = and/or to the axis Re{k)+Im{k) = 0. This additional 
symmetry is typical for lattices of size S x T, where both S and T are even numbers. In 
fact, it turns out that for this class of lattices the partition function Z{k) may be written in 
the following form: 

Z{k) = / n d^{l)d^{l) n F(jn) n crossin), (23) 

/ odd m odd " even 

where a given site (s, t) (with s = 1, . . . , S and t = 1, . . . , T) is said to be even or odd if the 
integer number s + 1 is respectively even or odd. While F{m) is the usual monomer term 
which we have introduced before in eq.(14), cross{n) is a new object obtained by multiplying 
the monomer term in the site n with the four dimer terms starting from or ending at the 
site n, and finally integrating with respect to the fermion fields in n: 

cross{n) = f d^jj{n)d^jj{n)F{n)Li{n)L2{n)Li{n - i)L2{n - 2). (24) 

By explicitely evaluating this expression, one finds that cross{n) may be written as 1 + 
k'^a{iljilj), where alipip) is a sum of products of four fermion fields defined in the neighbouring 
sites of n {i.e., n+ 1, n— 1, n + 2 and n — 2). Therefore the partition function Z{k) will be a 
polynomial in A;*^, with real (and positive) coefficients. As a consequence, if A; is a solution of 
Z{k) = 0, also ik will be a solution. In other words, the distribution of zeroes of Z{k) in the 
complex A;-plane will be invariant under the transformation {x,y) — > {—y,x): we will call 
this an I-symmetry. After combining this I-symmetry with the other P and C-symmetries, 
one immediately recognizes that the distribution of zeroes is invariant under refiections not 
only with respect to the real or the imaginary A;-axes, but also with respect to the axes 
Re{k) — Im{k) = and Re{k) + Im{k) = 0: this is in fact what one observes in [Fig.l]. 

With increasing lattice size the zero with the smallest imaginary part comes closer to 
the real axis. This zero of the partition function has been plotted in [Fig. 3] for various 
lattices of size S x T. As can be seen, for increasing T and fixed S the imaginary part of 
the zeroes approaches a finite, non-zero value, indicating that there is no phase transition 
on infinite stripes of extent S* < oo. We have then used a polynomial fit to derive, from the 
two previous series of points, the asymptotic values of the roots nearest to the real /c-axis 
corresponding to the lattices 2 x oo and 3 x oo. In practice, we have performed separate 
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fits of real and imaginariy parts of the sequence of roots for given S, using the polynomial 
ansatz, P{T) = ao + cii/T + a2/T^ + . . .. The number of coefficients in the polynomials was 
chosen so as to minimize the chi-squared value. It turns out that for the lattices with S = 2 
good fits are obtained if one uses a polynomial with six coefficients for the real part of the 
roots and a polynomial with five coefficients for the imaginary part. Instead, for the lattices 
with 5 = 3 good fits are obtained using polynomials with only three coefficients for both the 
real and the imaginary part of the roots. 

In this way we have found that the root nearest to the real fc-axis is approximately 
(0.464±0.003, 0.172±0.004) for the 2 x oo lattice and approximately (0.43±0.01, 0.09±0.01) 
for the 3 X cxD lattice. 

If one tries to linearly extrapolate these two points, assuming the existence of a real 
critical value kc in the thermodynamic limit S —>■ oo and T —>■ oo (see also below), one finds 
approximately (within the accuracy of our linear extrapolation) : 

fcc~0.39 . (25) 

Although the spatial extent of our lattices is still quite small, it is interesting to analyze the 
complex zeroes of infinite stripes under the assumption of a second order phase transition in 
the thermodynamical limit S — > oo. The real and the imaginary parts of the zero closest to 
the real axis are then expected to scale like: 

Re{k{S)) = kc+ " 



SI/U-' 



Im{k{S)) = -^. (26) 

Using this ansatz for our S = 2 and S = 3 results we find from the ^-dependence of 
Im{k) the critical exponent u ~ 0.63 and using this to extract kc from Re{k) we again find 
kc — 0.39. Of course, in order to get a realistic estimate of the errors on these numbers, one 
needs to consider also lattices with much larger S. 

This is surely an extremely interesting result: the value we have found for kc is consid- 
erably smaller than the value k = 1/2 we would have expected from naive argumentations. 
Since for j3 = oo (theory with free fermions, since e^ = 1/(3 -^ 0) it is well known that 
kc = 1/4 (in general it will be kc = l/2d, with d space-time dimensions), we are led by the 
result (25) to believe in the existence of a line of phase transition from (/? = 0,kc — 0.39) to 
(/? = oo, kc = 1/4). (See also ref. [Q, where this same result was found making use of the 
so-called eight-vertex model [^ for the analytical computations.) 

In order to get some further information about the question of the order of the phase 
transition, we have studied the behaviour of the (re-scaled) chiral condensate < ipip >, the 
(re-scaled) chiral susceptibility x ^^^ the (re-scaled) third derivative, with respect to 1/2/c, 
of the logarithm of the partition function Z in the vicinity of the critical point kc- The 
re-scaled chiral condensate < ijji(j > is given by the following expression in terms of the 
partition function Z{k) (see eqs.(5), (9) and (11)): 

7 7 1 ^ , /^N ^ k f 1 dZ{k)\ 



where A^ = S* x T is, as usual, the number of sites of the lattice. In [Fig. 4] we report the 
behaviour oi Q = — < ipip > as a function of k, for the two lattices 2 x 32 and 3 x 16. The 
re-scaled chiral susceptibility x is defined as: 



X 



(f 



MZ) 



{2kfNd(l/2kY 



-2 + 



k 

N 




{21 



In [Fig. 5] we report the behaviour of % as a function of k, for the same two lattices. It is 
easy to verify that x = —k{dQ/dk) — Q, so that x behaves as x — ~Q ~^ ~2 for A; ^ 0. 
Finally, we also consider the third derivative, with respect to 1/2/c, of the logarithm of the 
partition function Z: 



1 



d' 



V 



MZ) 



{2kfNd{l/2kf 
^ < ^jj{m)tjj{m)ilj{n)^jj{n)^jj{0)^jj{0) > -3N < ^jj^jj > x 



N^ < i)i) >' 



m,n 




(29) 



The behaviour of (^ as a function of k is shown in [Fig. 6], for the two lattices 2 x 32 and 
3 X 16. It is easy to demonstrate that (p = k{dx/dk) + 2%: from this one can immediately 
derive that ip behaves as (^ ~ 2% — ^ —4 for A; ^ 0. 

Clearly the three quantities are closely related. The sharpening of the crossover in 
< ipip > is reflected in the rise of the peak in x with increasing S, and the narrowing of 
the peak in x is expressed in terms of the rapidly rising peaks with opposite signature in 
(f. Certainly the behaviour of these three quantities is consistent with that expected for a 
continous phase transition, i.e. a second- or third-order phase transition. 

4 Conclusions 

We have evaluated analytically, using algebraic methods, the partition function Z for the lat- 
tice Schwinger model, with Nf = 1 Wilson fermions, in the strong-coupling regime {(3 = 0). 
For a given lattice SxT, the partition function is of the form Z{k, S, T) = {2k)~'^'^ Z{k, S, T), 
where A^ = S* x T is the total number of lattice sites and Z{k, S, T) is a polynomial in k of 
order 0{2N). By studying the zeroes in the complex plane {Re{k),Im{k)) of the partition 
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function Z{k, S, T) for a large series of small lattices S x T, we have found evidence for the 
existence of a critical value for the hopping parameter, which in the thermodynamic limit, 
S,T —>■ oo lies on the real axis at about kc — 0.39. We are led by this result to believe in 
the existence of a line of phase transition from [jS = 0,kc — 0.39) to [jS = cxd, kc = 1/4). In 
order to determine the order of the transition, it is clearly important to study in more detail 
the density of the zeroes near k^. This requires larger values of 5*. We have analyzed the 
chiral condensate < iptp >, the chiral susceptibility and the third derivative, with respect to 
1/2/c, of the partition function, in order to get some insights into the question of the order of 
the phase transition. Even though the present analysis does not yet allow to draw a definite 
conclusion on the order of the transition we have found some indications that the phase 
transition might be third order or even second order. 
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FIGURE CAPTIONS 

Fig. 1. The distribution of zeroes of the partition function Z{k, 2, 32) in the complex plane 

{Re{k),Im{k)). 

Fig. 2. The distribution of zeroes of the partition function Z{k, 3, 16) in the complex plane 
{Re{k),Im{k)). 

Fig. 3. The zeroes of Z{k,S,T) closest to the real A;-axis for various lattices of size S x T: 
crosses refer to lattices with S = 2 and T = 2, 3, 4, 6, 8, 10, 16, 32, while triangles refer 
to lattices with 5 = 3 and T = 3, 4, 5, 6, 8, 9, 10, 12, 16. 

Fig. 4. The re-scaled chiral condensate, Q = — < ipip >, as a function of k, for the two lattices 
2 X 32 and 3 X 16. 

Fig. 5. The re-scaled chiral susceptibility x (defined in the text by eq.(28)) as a function of 
k, for the two lattices 2 x 32 and 3 x 16. 

Fig. 6. The re-scaled third derivative (p{k) of ln{Z) with respect to 1/2A; (defined in the text 
by eq.(29)), for the two lattices 2 x 32 and 3 x 16. 
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Table 1: The coefficients of the polynomial Z{k, 2, 32) = J2^=o a4„A;*^". 

n a^n (for the 2 x 32 lattice) 

1 

1 512 

2 129024 

3 21331968 

4 2602369024 

5 249774997504 

6 19638234644480 

7 1300255043747840 

8 73961277759684608 

9 3668969473236271104 

10 160569828529865228288 

11 6255905737448198504448 

12 218535073687128684625920 

13 6883374618372455350140928 

14 196354147234439285612478464 

15 5089687819224732967794376704 

16 120172957550769852891363540992 

17 2588628553768726317536821379072 

18 50910684319442948897317601673216 

19 914109924603780980467032270045184 

20 14970604113895231454059152342515712 

21 223213979431827944404264819095502848 

22 3021076338904294430787178751238078464 

23 36958171640461419289379070931541950464 

24 406259327271308991154552486860991496192 

25 3980443123196403924685803149026276671488 

26 34377853959664866223457481686909255155712 

27 257704138028993570924127408162720968605696 

28 1639689392947634033166698976649965853474816 

29 8559749426407522788060672587795634988253184 

30 34663799157267802700683573313175561579790336 

31 97874256443064108849693716776331356231696384 

32 146811384664566690398085268830912235998019584 



48 ^ u2n 



Table 2: The (non-zero) coefficients of the polynomial Z(k, 3, 16) = J2'n=o '^2nk 



n 


o-2n (for the 3 x 16 lattice) 



3 


1 
1024 


4 


768 
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6144 
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516096 
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811008 
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6334464 
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174620672 


10 


431947776 
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3203923968 


12 


45588021248 
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153911033856 


14 


1063856898048 


15 


9906042175488 
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40917798223872 
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261009935695872 
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1857033179496448 
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2038448818744602390429696 
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8604050273287356083601408 
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3492308812582291 1234703360 
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136172615876993812147470336 
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515227117640733139176259584 


36 


1857978729092258641800593408 


37 


6360979177485248012381847552 


38 


20656762024719021784826904576 


39 


62433587250829914862254030848 


40 


174406412859094327618009300992 


41 


446095973409032980536257150976 


42 


1025256020010205168904800567296 


43 


2088088022191744886698175102976 


44 


3738100459121645858804962689024 


45 


5610148109882451164480363560960 


46 


7107243167143739985776879861760 


47 


6755399440834383071 1 15485380608 
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